Interactions, disorder and local defects in graphite 
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Recent experiments report the existence of ferromagnetic and superconducting fluctuations in 
graphite at unexpectedly high temperatures. The interplay of disorder and interactions in a 2D 
graphene layer is shown to give rise to a rich phase diagram where strong coupling phases can 
become stable. Local defects can explain the ferromagnetic signals. 
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Introduction. A number of recent experiments suggest 
that pure graphite behaves as a highly correlated elec- 
tron system[l|. In particular it shows a metal-insulator 
transition in magnetic fields and insulating behavior in 
the direction perpendicular to the planes in different 
samples [2i,_33 A.:.y>-,3.,l,AiM- The interest in this material 
is focussed nowadays in the observation of ferromagnetic 
behavior ^3 J enhanced by proton bombardment 'l l] what 
opens a new way to the creation of organic magnets. 

In refs. [H m a simple microscopic model was pro- 
posed as a new framework to study the physics of 2D 
graphene sheets and its topological variant fuUerenes 
and carbon nanotubes. The main assumption of the 
model is to neglect the coupling between layers and 
consider graphite as a pure two-dimensional system. 
This assumption is supported by experiments where an 
anisotropy of up to three orders of magnitude is mea- 
sured in magnetotransportj^. The model predicts non- 
Fermi liquid behavior for the graphene system and can 
account for the linear behavior with energy of the quasi- 
particle scattering rate|l4| observed in photoemission 
experiments 0. 

In this work we review the main features of the model 
with and without disorder and propose a new mechanism 
to explain the ferromagnetic fluctuations observed in the 
experiments. 

The model. RG results. The conduction band of 
graphite is well described by tight binding models which 
include only the tt orbitals which are perpendicular to the 
graphite planes at each C atomT^. The two dimensional 
hexagonal lattice of a graphene plane has two atoms per 
unit cell. A tight binding calculation with only nearest 
neighbors hopping gives rise to the dispersion relation 
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£;(k) = ±ty 1 -I- 4cos2 —kx -I- 4 cos —kx cos -ky (1) 

whose lower branch is shown in Fig^ This dispersion 
relation gives rise at half filling to a Fermi surface con- 
sisting of six isolated points two of which are inequiva- 
lent. A low-energy effective Hamiltonian can be defined 
by expanding the dispersion relation about any of the 
Fermi points. The resulting Hamiltonian has the form 
of a massless two dimensional Dirac Hamiltonian. The 




FIG. 1: Lower branch of the electronic dispersion relation. 
The cusps appear at the six corners of the first Brillouin zone. 



Fermi velocity, vp, can be expressed in terms of the ma- 
trix elements between nearest neighbor tt orbitals, t, as 
Vp = (3ta)/2, where a is the C-C distance. 

The electronic states within each graphene plane are 
described by two two-component spinors associated to 
the two inequivalent Fermi points in the Brillouin Zone. 
The Hamiltonian of the free system is: 



Hq = ivp / d^j:\E'(x)7 • V*(f) 



(2) 



where the two-dimensional 7 matrices are built as appro- 
priate combinations of Pauli matrices 12i. The Hamilto- 
nian Q gives an effective description of graphite in an 
energy range bound by a lower cutoff .27eV dictated 
by the interlayer coupling, and a higher cutoff, where the 
bands can no longer be approximated by a linear disper- 
sion relation ^ 3 — 4eV. 

The Hamiltonian (O is the perfect model for Renor- 
malization Group (RG) calculations. It is scale invari- 
ant and does not have the complications of an extended 
Fermi surface. The model is similar to the D — 1 electron 
system in that it has Fermi points and linear dispersion. 
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Nevertheless naive dimensional analysis shows that four 
and more Fermi interactions are irrelevant in this two- 
dimensional case. The only interaction that may survive 
at low energies is the long (infinite) range Coulomb in- 
teraction, unscreened because of the vanishing density 
of states at the Fermi point. Following the quantum 
field theory nature of the model, we trade the classical 
Coulomb interaction 

(3) 

where g = jAi^vp is the dimensionless coupling con- 
stant, by a local gauge interaction through a minimal 
coupling. 

Lint =9 J d?xdtf{x, t)A^,{x, t) , (4) 
where the electron current is defined as 

This interaction is marginal in the RG sense, all the rest 
are irrelevant. The RG analysis of the model gives the 
following results [T^IT^: 

1. From the computation of the electron self-energy at 
the one loop level we get a non trivial renormalization of 
the Fermi velocity that grows in the infrared. This 
result implies a breakdown of the relation between the 
energy and momentum scaling, a signature of a quantum 
critical point. 

2. The electron-photon vertex and the photon prop- 
agator are not renormalized at the one loop level. This 
means that the electric charge is not renormalized, a re- 
sult that could be predicted by gauge invariance, and 
it also implies that the effective coupling constant g — 
e^/iirvp decreases at low energies defining an in- 
frared free fixed point of the RG. It is interesting to note 
that the Lorentz invariance of the model that was ex- 
plicitly broken by the Fermi velocity is recovered at the 
fixed point since the velocity of light, c, fixes a limit to 
the growing of the Fermi velocity. 

2. From the electron self-energy at two loops order we 
get a non trivial wave function renormalization mean- 
ing that the infrared stable fixed point corresponds to 
a free fixed point different from the Fermi liquid. This 
result has been shown to persist in the non-perturbative 
regime ^3- This is a non-trivial result that has physi- 
cal implications. In particular it implies that the inverse 
quasiparticle lifetime increases linearly with ener gy [I^ . 
a result that has been observed experimentally in '15*1 in 
the energy range of validity of the model. 

In conclusion, we have shown that without disorder, 
edges, or other perturbations, the graphene system at low 
energies has gapless excitations differing from the Fermi 
liquid quasiparticles but does not support magnetic or 
superconducting instabilities. 




FIG. 2: Formation of a pentagonal ring in the honeycomb 
lattice. Points a,b,c,d... have to be identified with points 
a', b' , c' ,d' .... The defect can be seen as a disclination, defined 
by the straight dashed lines. 



The strong coupling regime of the graphene system 
has been analyzed in Il3 ■ There it is argued that a 
dynamical breakdown of the chiral symmetry (degener- 
acy between the two Fermi points) will occur at strong 
coupling and a gap will open in the spectrum forming a 
kind of charge density wave. Graphite can then be seen 
as an excitonic insulator that can become ferromagnetic 
upon doping. Being non-perturbative, these phenomena 
are likely to be washed out by any amount of disorder at 
intermediate energies. 

Inclusion of disorder. 

The previous description analyzes the small momen- 
tum scattering due to the long range Coulomb interac- 
tion, as it is the only one which leads to logarithmically 
divergent perturbative corrections. Some electronic in- 
stabilities, like ferromagnetism or anisotropic supercon- 
ductivity, require the existence of short range interac- 
tions with significant strength. The irrelevant character 
of short range interactions can be changed by the pres- 
ence of disorder that enhances the density of states at 
the Fermi level. 

Disorder can be included in the renormalization group 
scheme by the introduction of random gauge fields. This 
is a standard procedure in the study of the states de- 
scribed by the two dimensional Dirac equation asso- 
ciated to random lattices or to integer quantum Hall 
transitions "20*, '21', '23|. There it is seen that, usually, the 
density of states at low energies is increased. To demon- 
strate how these gauge fields can arise in the graphene 
system, we will describe in detail a special type of disor- 
der that we call topological disorder. 

The formation of pentagons and heptagons in the lat- 
tice, without affecting the threefold coordination of the 
carbon atoms, lead to the warping of the graphene sheets, 
and are responsible for the formation of curved fuUerenes, 
like Geo. They can be viewed as disclinations in the lat- 
tice, and, when circling one such defect, the two sub- 
lattices in the honeycomb structure are exchanged (see 
Fig- El)- The two fermion flavors associated to the two 



Fermi points are also exchanged when moving around 
such a defect. The scheme to incorporate this change in 
a continuum description was discussed in^3|- The pro- 
cess can be described by means of a non Abelian gauge 
field, which rotates the spinors in flavor space. The vec- 
tor potential is that of a vortex at the position of the 
defect, and the flux is ±7r/2. 

Dislocations can be analyzed in terms of bound discli- 
nations, that is, a pentagon and an heptagon located at 
short distances, which define the Burgers vector of the 
dislocation. Thus, the effect of a dislocation on the elec- 
tronic levels of a graphene sheet is analogous to that of 
the vector potential arising from a vortex-antivortex pair. 
We can extend this description 23] , and assume that a 
lattice distortion which rotates the lattice axis can be 
parametrized by the angle of rotation, 6{r), of the local 
axes with respect to a fixed reference frame. Then, this 
distortion induces a gauge field such that: 

A(r)=3V^(r)(^^ ) (5) 

Thus, a random distribution of topological defects can be 
described by a (non abelian) random gauge field. 

Other types of disorder can similarly be associated to 
random gauge fields. The complete Hamiltonian of the 
system can be written as 

H = Hee + Hdisorder (6) 

where 

(7) 

Hd^sorder j {x)T^ {x)A{x) (8) 

vr characterizes the strength and the 4x4 matrix T 
the type of the vertex. In general, A{x) is a quenched, 
Gaussian variable with the dimensionless variance A, i.e., 

{A{x))=0 , {A{x)A{x')) ^ AS^{x-x') . (9) 

In ref. (2^ a complete RG study of the disordered system 
was analyzed by adding gauge couplings associated to 
all possible gamma matrices, i) For a random chemical 
potential, the 4x4 matrix F is given by F = 79. The 
long range components of this type of disorder do not 
induce transitions between the two inequivalent Fermi 
points. This type of disorder yields an unstable fixed 
line, ii) A random gauge potential involves the 4x4 
matrices F — iji and F = 172. This type of disorder gives 
rise to a stable fixed line which is linear in the {g, A)- 
plane. iii) (a) A fluctuating mass term is described by 
F = 14x4. (b) Topological disorder is given by F = 175 
with 75 = l2x2 ^ o'2. c) To complete the discussion, we 
also mention F = 175 where 75 = l2x2 (Si cri- This vertex 
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FIG. 3: One-loop phase diagram for two-dimensional mass- 
less Dirac spinors including long-ranged electron-electron in- 
teraction g and disorder A. Top: Random chemical potential 
(r — 70). Center: Random gauge potential (F = 171,172). 
Bottom: Random mass term (F = 14x4), topological disor- 
der (r = 175), and r — 175. 



type can be related to an imaginary mass that couples 
the two inequivalent Fermi points. All these types of 
disorder will yield a stable fixed line which is cubic in 
the {g, A)-plane. 

The phase diagram obtained in |24| | is reproduced in 
Fig. 3. 

i) For a random chemical potential (F = 70), wr = 
vi remains constant under renormalization group trans- 
formation. There is an unstable fixed line at Vp = 
vfA/ {2e'^). In the (g. A) -plane, the strong-coupling and 
the weak-coupling phases are separated by a hyperbola, 
with the critical electron interaction g* = /v*p — 
2e*/(w^A). ii) A random gauge potential involves the 
vertices F = 171 , 272 . The vertex strength renormal- 
izes as = vp. There is thus an attractive Luttinger- 
like fixed point for each disorder correlation strength A 
given by v*p = 2e^/A or g* — A/2, iii) For a random 
mass term F = l4x4, topological disorder F = 175, and 
r = *75, we have vr = Vp/v^. There is thus again 
an attractive Luttinger-like fixed point for each disor- 
der correlation strength A given by Vp — y^2v^e^/A or 
g* = i/Aey{2vi). 

The most interesting phase is the one induced by a ran- 
dom gauge potential, a random mass term or the topo- 
logical disorder. All of them drive the system towards 
a new stable, Luttinger-like fixed point. This phase is 
characterized by a vanishing quasiparticle residue, lead- 
ing to anomalous one particle properties. The Luttinger 
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liquid features associated to this fixed line are notoriously 
difficult to observe, although they can be probed in tun- 
neling experiments, or by measuring the peak width in 
ARPES. They will also influence the interlayer transport 
properties 25, 26] . Small perturbations by other types 
of disorder, like a random local potential induce a flow 
along this fixed line, as in the absence of interactions '2^. 
The strong coupling fixed point describes, most likely, a 
disordered insulating system. 

Localized states. In addition to the extended disorder 
discussed previously, a gra.phcnc plane can show states lo- 
calized at interfaces [28l l29| . which, in the absence of other 
types of disorder, lie at the Fermi energy. Changes in the 
local coordination can also lead to localized states 30J. 

The tight binding model defined by the tt orbitals at 
the lattice sites can have edge states when the sites at the 
edg e belong all to the same sublattice(zig-zag edges') psl 
l29l |. These states lie at zero energy, which, for neutral 
graphene planes, correspond to the Fermi energy. 

In a strongly disordered sample, large defects made up 
of many vacancies can exist. These defects give rise to 
localized states, when the termination at the edges is lo- 
cally similar to the surfaces discussed above 31]. Note 
that, if the bonds at the edges are saturated by bond- 
ing to other elements, like hydrogen, the states at these 
sites are removed from the Fermi energy, but a similar 
boundary problem arises for the remaining tt orbitals. A 
particular simple example is given by the crack shown in 
Fig.m. 

These states are half filled in a neutral graphene plane. 
In the absence of electron interactions, this leads to a 
large degeneracy in the ground state. A finite local re- 
pulsion will tend to induce a ferromagnetic alignment of 
the electrons occupyin g th ese states, as in similar cases 
with degenerate bands pj3|. Hence, we can assume that 
the presence of these states leads to magnetic moments 
localized near the defects. 




FIG. 4: Example of a crack in a graphene plane. The atoms at 
the upper edge and those at the lower edge belong to different 
sublattices. 



We now have to analyze the influence of these mag- 
netic moments in conduction band described in the pre- 
vious sections. The hopping between the states involved 
in the formation of these moments and the delocalized 
states in the conduction band vanishes by definition, if 
the localized states lie at zero energy. Hence, a Kondo 
like coupling mediated by the hopping will not be in- 
duced. The localized and conduction band states, on the 
other hand, are defined on the same lattice sites. The 
existence of a finite local repulsion, U, will lead to an 
effective ferromagnetic coupling. The local coupling, at 
site i, between the localized states and the conduction 
band is proportional to U J2j Pi,j ; where pij is charge of 
state j at site i. In order to get an order of magnitude 
estimate of the effect of these states, we will assume that 
the number of states induced near a vacancy is similar to 
the number of atoms at its edge, Af, and that these states 
are sufficiently localized around the vacancy. Hence, each 
vacancy nucleates a moment of order Af. The effective 
coupling between a vacancy and the conduction electrons 
is proportional to UAf^ and it is distributed over an area 
^ A/"^ . The conduction electrons will mediate an RKKY 
interaction between the localized moments: 

jRKKvir) ^U' f d^ke^i^^xlk) ~ U^^j^ (10) 
J VF\r\-^ 

Where the static susceptibility is xi^) l^l'l^, and a 
is the lattice constant. It is interesting to note that, due 
to the absence of a finite Fermi surface, the RKKY in- 
teraction in ea. H10(l does not have oscillations. Hence, 
there are no competing ferro- and antiferromagnetic cou- 
plings, and the magnetic moments will tend to be fer- 
romagnetically aligned, leading to an effective magnetic 
field, Hext(j): with non zero average, acting on the con- 
ducting electrons. 

From power counting, this coupling is relevant in the 
Renormalization Group sense. Thus, in the presence of 
extended vacancies, the RG flow discussed in the previ- 
ous sections has to be arrested at scales comparable to 
{Hext{Y)) ^ UAfpvac, whcrc pvac is the concentration of 
the large vacancies which may give rise to localized states. 
At lower energies, or temperatures, the graphene planes 
with extended vacancies will behave as a ferromagnet. 

Conclusions. In this work we present a microscopic 
model for studying the low energy properties of a sin- 
gle graphene layer as a model relevant for some graphite 
samples showing two-dimensional anomalous behavior. 
In particular we tried to envisage a model able to ex- 
plain the ferromagnetism observed recently in a variety 
of graphitic materials. 

The model is based on the particular dispersion re- 
lation of the 2D honeycomb lattice that, at half filling, 
has Fermi points instead of Fermi lines. The lineariza- 
tion of the dispersion about a Fermi point gives rise to 
a model similar to the one-dimensional electron system 
with zero density of states at the Fermi level. Unlike the 
ID case, the real two-dimensional nature of the present 
model makes the four Fermi interactions irrelevant in 
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the rcnormalization group sense while the long range 
Coulomb interaction is unscreened and plays an impor- 
tant role. It renormalizes the Fermi velocity that grows 
at low energies while the effective charge is not renormal- 
ized, a consequence of the gauge invariance. The effective 
coupling constant jvp goes to zero driving the system 
to a non-trivial infrared free fixed point. As a conse- 
quence of the singular Coulomb interaction, the electron 
acquires anomalous dimension and the quasiparticle scat- 
tering rate grows linearly with frequency at intermediate 
frequencies, as observed in experiments. The model does 
not support magnetic or any other short range interac- 
tions at this level. 

The presence of disorder changes the previous situa- 
tion in various respects. We have considered two types 



of disorder, non-local disorder as the one produced by im- 
purities or lattice distortions, modelled by the coupling of 
the electrons to random gauge fields as in the random lat- 
tice models, and local large defects as the ones produced 
in the experiments by proton bombardment. Extended 
disorder gives rise to a rich phase diagram with strong 
coupled phases whose physical properties are still to be 
analyzed. Local defects give rise to the appearance of 
local moments whose interaction can induce ferromag- 
netism in large portion of the sample. 
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acknowledged. 
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